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SIMILARITY  PARAMETERS  FOR  RADIATIVE  ENERGY  TRANSFER 
IN  ISOTHERMAL  AND  NON -ISOTHERMAL  GAS  MIXTURES* 

by 

S.  S.  Penner,*  M.  Thomas**  and  G.  Adomeit*** 

Karman  Laboratory  of  Fluid  Mechanics  and  Jet  Propulsion 
Daniel  and  Florence  Guggenheim  Jet  Propulsion  Center 
California  Institute  of  Technology 
Pasadena,  California 

ABSTRACT 

The  similarity  groups  for  multicomponent,  reacting 
gas  mixtures  with  radiative  energy  transport  are  derived 
(Section  I).  The  resulting  relations  are  used  to  consider 
the  feasibility  of  scaling  for  flow  processes  with  radiative 
energy  transport  under  highly  simplified  conditions 
(Sections  II  and  III).  Next  the  scaling  parameters  are 
derived  for  radiant  energy  emission  from  isobaric  and 
isothermal  gases  for  arbitrary  opacities  and  various 
spectral  line  and  molecular  band  models  (Section  IV). 

Scaling  parameters  for  radiant  energy  emission  from 
isobaric  but  non-isothermal  systems  are  discussed  for 
arbitrary  opacities  and  various  spectral  line  and  molecular 
band  models  under  the  restrictions  imposed  on  the  allowed 
temperature  profiles  for  dispersion  and  Doppler  lines  by 
the  Eddington -Bar bier  approximation  (Section  V).  Finally, 
we  consider  the  radiative  scaling  properties  for  represen¬ 
tative  temperature  profiles  for  both  collision-broadened 
and  Doppler-broadened  line  profiles  on  the  basis  of  exact 
numerical  calculations  that  we  have  performed  for  a 
rotational  spectral  line  belonging  to  a  molecular  vibration- 
rotation  band  (Section  VI).  It  appears  that  simple  scaling 
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rulea  generally  constitute  a  fair  approximation  for 
dispersion  lines  in  non-isothermal  systems  but  that 
corresponding  relations  apply  to  lines  with  Doppler 
contour  only  in  the  transparent  gas  regime. 


I.  DETERMINATION  OF  SIMILARITY  PARAMETERS 

The  techniques  for  identifying  the  similarity  groups  for  systems 
described  by  a  set  of  conservation  equations  are  well  known  and  have 
been  developed  in  detail  previously  for  reacting,  multicomponent  gas 
mixtures.  ^  Although  the  particular  form  of  the  conservation  equations 
used  in  the  original  analysis  has  been  shown  to  be  in  error  in  several 
respects,2'^  the  resulting  similarity  groups  are  unchanged  if  the  same 
approximations  are  used  as  in  the  original  article.  *  Some  generali¬ 
zation  of  procedure  can  be  made  without  difficulty  by  following  the  method 
described  below  j^Eqs.  (4)  and  (6)  are  not  given  in  ref.  lj  . 

We  start  with  the  set  of  conservation  equations  jjsee  Eqs.  (12), 
(13),  (16)  to  (18),  and  (23)  in  Chapter  XVIII  of  lj  ;  next  we  use  a  binary 
mixture  approximation  for  each  of  the  diffusion  velocities  fusing  Eq.  (4) 
on  p.  243  of  I,  replace  the  subscript  1  by  K  (K=l,  2,  .  .  .  n)  for  an  n- 

_l 

component  mixture,  replace  by  1-YK  and  W^  by  W  =  average  molecular 

*  S.S.  Penner,  Combustion  Researches  and  Reviews  1955,  pp-  140-162, 
Butterworths  Scientific  Publications.  London  195b . 

2  For  a  derivation  of  the  relevant  equations  from  continuum  theory,  see 
W.  Nachbar,  F.A.  Williams  and  S.S.  Penner,  Quart.  Appl.  Math.  17_, 

43  (1959). 

The  correct  results  are  also  given  by  S.S.  Penner,  Chemistry 
Problems  in  Jet  Propulsion,  Pergamon  Press,  Ltd.,  London  1957. 

This  book  will  be  referred  to  as  I  hereafter. 
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w eight  of  the  fluid  mixture  remaining  without  species  Kj  ;  finally  we 
write  the  complete  heat  flux  vector  j^see  Eq.  (23)  on  p.  239  of  ij  in  the 
form^ 

«  ■  •  [»♦  >ch] 

Here  A  is  the  thermal  conductivity  associated  with  molecular  collisions 

and  A  .  identifies  the  thermal  conductivity  associated  with  chemical 
ch 

4 

reaction;  A  .is  given  explicitly  by  the  expression 
cn 

\h  -t  t  DKK'WK'HK  TT 

K=1  K'=l 

where  (c)  =  total  number  of  moles  per  unit  volume  of  mixture,  f>  s  fluid 
density,  DRKi  =  multicomponent  diffusion  coefficient  which  is  a  known 
function  of  the  binary  diffusion  coefficients  and  the  mixture  composition,  5 
WK,=  molecular  weight  of  species  K',  HK  =  molar  enthalpy  of  species  K, 
X^,  =  mole  fraction  of  species  K\  T  =  temperature. 

Using  the  specified  starting  relations  and  standard  procedures,1 
we  find  the  following  set  of  similarity  parameters  for  multicomponent, 
reacting  gas  mixtures  without  radiative  energy  transport  if  the  subscript 
o  identifies  suitably  chosen  reference  conditions; 

Yq  a  ratio  of  specific  heat  at  constant  pressure  (  cp  Q)  to  the  (1) 

specific  heat  at  constant  volume  ( c  )  for  the  fluid 
r  V,  o 

mixture, 

J.O.  Hirschfelder ,  University  of  Wisconsin  Report  WIS-ONR-18, 
February  6,  1956;  J.N.  Butler  andR.S.  Brokaw,  J.  Chem.  Phys.  26, 
1636  (1957). 

5  J.O.  Hirschfelder,  C.F.  Curtiss  and  R .  B.  Bird,  Molecular  Theory 
of  Gases  and  Liquids,  John  Wiley  and  Sons,  Inc. ,  New  York  19^4 . 
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Reynolds  numbers  a  Pe.  =  ^v^L. //iQ  (f  =  density,  v  =  flow  (2) 

velocity,  L.  =  ith  characteristic  length,  jt.  =  mixture 
viscosity;  i  =  1,  2,  . .  .  J  ), 

Schmidt  numbers  a  Sc..  -  ji  /  o  D  ( D  =  diffusion  (3) 

IN.  O  'O  IN,  O  K,  O 

coefficients  for  species  K;  K  =  1,  2,  . .  .n). 


Sc  1  3  >u  / p  k  D  (k  =  thermal  diffusion  ratio  for  (4) 

K  o  ) o  IK,o  K,o  TK 

species  K,  D-..  =  k__,  D  =  thermal  diffusion 

1  IN. »  O  I  IN,  O  ,  O 

coefficient  for  species  K), 


Prandtl  number  a  Pr  =  c  M  /  A  , 

p,  o  o  o 

Pr*  =  c  Ml  A  ,  (this  group  will  actually  remain  invariant 
p, o  o  ch, o  ®  r  7 

if  the  groups  Pr,  Sc  and  D  .  are  fixed), 
in  111 ,  i ,  r 


Mach  number  5  M  =  y  f0y0  /r0PQ  ^  =  Pre88urc)» 

Froude  numbers  as  v^/gL.  (g  =  gravitational  acceleration), 

Damkohler's  first  similarity  groups  =  D,  .  =  L.U  /v 

76  r  I,i,  r  1  r,  o  o 

(Uf  ^  =  characteristic  reaction  frequency  for  the  rth 
chemical  process;  r  =  1,  2,  . .  .  m). 


(5) 

(6) 

(7) 

(8) 
(9) 


Damkohler's  third  similarity  groups  2  DITI  .  =  q  'U  L. /v  c  T 

76  r  III,  i,  r  nr  r ,  o  l  o  p,o  o 

(q  '  =  heat  release  per  unit  mass  in  the  rth  chemical  reaction), 

r  (10) 


4  =  i  v  2/(c  /r  ) t 

J  2  o  p,  o  o'  c 


(ID 


A.  Radiative  Energy  Transport  in  the  Diffusion  Approximation 

Throughout  the  following  discussion  we  neglect  photochemical 


For  an  elementary  discussion  of  the  basic  equations,  see  S.S.  Penner 
and  R.W.  Patch,  "Radiative  Transfer  Studies  and  Opacity  Calculations 
for  Heated  Gases",  Technical  Report  No.  6,  Contract  AF  49(638)-984, 
California  Institute  of  Technology,  Pasadena,  Calif.  ,  January  1962;  in 
press  in  Papers  Presented  in  Honor  of  Modesto  Panetti,  Politecnico  di 
Torino,  Accademia  Nazionale  dei  Lincei,  Rome,  Italy,  1962. 
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reactiona.  In  the  diffusion  approximation,  the  effective  thermal 
conductivity  ia  augmented  by  the  term  Ara  where 


Ari  =  .6rT3/3kLiBo. 

V  =  Stefan-Boltzmann  constant  and  =  Roaaeland  mean  absorption 

coefficient.  Hence  an  additional  analogue  of  the  Prandtl  number  will 
appear,  viz. , 


Pr*  =  c  Ml  A 

p.o  o  ra,o 


(12) 


In  the  momentum  equation,  an  additional  body  force  per  unit 


volume 

r  a  =  .(16kT5/3c)7T 

appears,  where  c  =  velocity  of  light.  This  body  force  per  unit  volume 
may  be  considered  as  an  additive  term  to  the  pressure  gradient,  i.e., 
the  fluid  pressure  is  replaced  by  the  sum  of  the  fluid  and  radiation 
pressures.  Thus  a  new  similarity  group  may  be  formed  which 
measures  essentially  the  ratio  of  these  pressures.  It  may  be  written  as 


16  G"T  4  ,  ,  (4  •'T  4  /  c ) 

o  Z  Y  tZ  o 

5cprt  3  o  7l  2 1 


(rfovo  > 


The  ratio  of  the  radiant  energy  density  for  a  biackbody  at  the  reference 

4 

temperature  Tq  (=4^Tq  /c)  to  the  translational  energy  density  under 
reference  conditions  (=/0v0^/2)  designate  as  the  similarity  group 


4rT  4/c) 

«  s  °. 

i  rz 

p  v  Z 
Too 


(13) 


The  preceding  considerations  indicate  that,  in  the  diffusion 
approximations,  allowance  for  radiative  energy  transfer  leads  to  the 
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requirement  that  an  additional  analogue  to  the  Prandtl  number  ^  see 

Eq.  (12)j  and  the  aimilarity  group  ^  be  considered.  The  reciprocal 

of  the  Rosseland  mean  absorption  coefficient  is  the  Rosseland  mean 

free  path  1_  and  is  defined  by  the  relation 
r  Ro  ' 


^Ro  =  ^L,Ro) 


-i.i5  r 


i 


4  x 
x  e 


L,V,T  (e  -1) 


dx, 


hv 
x  =  P f 


where  k  identifies  the  linear  absorption  coefficient  in  the 

V  t  I 

frequency  interval  between  V  and  v  +  dv  at  the  temperature  T.  Thus 
constancy  of  the  similarity  group  Pr*  imposes  severe  restrictions  on 
the  allowed  values  of  k.  ^  _  although  it  is  not  required  that  the  products 
of  the  linear  spectral  absorption  coefficients  and  characteristic  lengths 
be  maintained  invariant,  as  is  necessary  for  radiative  transfer  problems 
for  arbitrary  opacities  (see  Section  IC)  in  isothermal  systems. 


B.  Radiative  Energy  Transport  for  Transparent  Gases 

The  radiation  pressure  for  transparent  gases  without  external 
source  is  negligibly  small.  For  isothermal  systems,  the  energy  loss 
from  unit  volume  by  radiation  is^ 


7.  rim*«rT4ELiPI 

where  E.  is  the  Planck  mean  absorption  coefficient  which  equals 
L i  Jr*l 

*♦ 

the  spectral  emissivity  per  unit  length.  The  term  7 .  Fr^  occurs  as 
an  additional  term  to  the  internal  energy  flow  rate  per  unit  volume 

^  v  .  7  u 

where  u  represents  the  specific  internal  energy.  Hence  we  may 


construct  the  ratio 
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‘4<rT4t‘  /rT4t.  D1  L.  \r(T*)4t.  * 

_ L,P1  _it  O  L.Pl.o  i  L,  PI 


p  v  .  vu 


P  v  c  T 
f o  o  p,  o  o 


*  _*  *  4 

f  V  .  Vu 


where  the  starred  quantities  are  dimensionless.  The  parameters 

<TT  4t.  L.  - 

r  = _ (i 

1  P  v  c  T 

JO  o  p,  o  o 

measure  the  ratio  of  radiative  energy  loss  from  the  system  per  unit 
surface  area  to  the  free  stream  rate  of  enthalpy  transport  per  unit 
area.  The  parameter  T  .  plays  an  important  role  in  problems  on 
stellar  turbulence.^ 

Other  important  similarity  groups  in  flows  with  radiant 
energy  transfer  are  the  Bouguer  numbers 

BuV,i"  kL,V,T,oLi  (1 


BuPl,i  =  ^L,Pl,oLi 


(15a) 


Bu_  .  *5  L. 
Ro,i  L,Ro  i 


(15b) 


and  the  Boltzmann  number 


—  -  Q 

r  v  c 
'  o  o  p, o 

The  Boltzmann  number  is  seen  to  arise  from  the  similarity  group  P  . 
by  dividing  this  quantity  by  the  ith  Bouguer  number  BUpj  .. 


A.  Unsold,  Physik  der  Stcrnatmospharcn,  Julius  Springer,  Berlin  1958 
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C.  Radiative  Energy  Transport  for  Arbitrary  Opacitiea* 

When  neither  the  diffusion  approximation  nor  the  transparent 

gas  approximation  are  applicable,  proper  allowance  for  radiative  energy 

transport  becomes  particularly  difficult  since  local  radiant  energy 

contributions  are  determined  by  integrals  over  the  accessible  field  of 

view.  Neglecting  photochemical  processes,  it  appears  now  that 

complete  simulation  is  possible,  even  for  isothermal  systems,  only 

if  a  set  of  similarity  groups  involving  the  parameters  Bu  .  remains 

v,  1 

invariant  for  all  important  lengths  L.  at  all  frequencies  V  .  For  non- 
isothermal  systems,  the  problem  becomes  complicated  still  further 
because  simulation  now  requires  invariance  of  similarity  groups 
involving  the  spectral  radiant  flux  per  unit  solid  angle  at  x^  with  the 
solid  angle  measured  (see  Fig.  1)  in  the  direction  (3^ -3j)/ |  x., **i  • 
i. e. ,  along  C,  viz.  , 


where  the  symbol  C  indicates  that  integration  is  to  be  performed 
along  the  path  defined  by  a  straight  line  drawn  between  the  ends  of 
(initially  chosen)  vectors  x.,  and  Xj.  Here  the  scalar  spectral  black- 
body  radiancy  R°  Q  and  the  scalar  linear  spectral  absorption  coefficients 
are,  in  general,  complicated  functions  of  the  special  location. 

This  important  problem  is  considered  in  detail  in  Section  IV  for  iso¬ 
thermal  emitters  and  in  Section  V  for  non -isothermal  systems. 
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Fig.  1.  Schematic  diagram  showing  the  definition*  of 
x,  and  C 


II.  SCALING  OF  FLOW  WITH  MAINTENANCE  OF  CONSTANT  VALUES 
OF  D.  ,  AND  r.  FOR  TRANSPARENT  GASES 

An  interesting  special  case,  which  is  of  some  practical 
importance,  involves  the  scaling  of  rsactive  gas  flows  for  a  given 
chemical  system  in  such  a  way  that  the  reference  temperature,  velocity 
and  specific  heats  are  invariant. 

Let  us  identify  by  the  subscripts  H  and  Mo  a  large-scale  burner 
and  the  model,  respectively.  It  follows  then  from  Eqs.  (9)  and  (14)  that 


Ll;HUr,  o;H  _  Li;MoUr,  o;Mo 


o;H 


ojMo 


and 


<rTo;H3^L,Pl;o;HLi;H  _fTo;MQ3  *L,  PI,  o;MoLl;Mo 


f  o;Hvo;Hcp,  o;H 


f  o;Movo;Mo  p,  o;Mo 
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"n  a  5f~  "7> 

'o:Mo 

denote  the  scaling  parameter  for  the  physical  variable  ^  .  Furthermore. 

let  us  consider  two  propellant  streams  that  are  injected  at  the  same 

reference  conditions  with 

n  =  n .  =  n_  =  1. 
v  c  T 
P 

The  specified  requirements  for  similarity  now  reduce  to  the  relations 


Since 


nL  nU  =  1 
i  r 


i  n  =  n,  . 

E  Li  / 

L,P1 


nu  *  "f 


if  Qf  represents  the  overall  order  of  the  rth  chemical  process,  we  may 
rewrite  Eq.  (18)  in  the  form 


V  nf 


Therefore,  the  similarity  groups  D.  .  will  remain  invariant  if  model 

tests  are  performed  on  small-scale  engines  at  elevated  density  since  &r 

will  generally  be  larger  than  unity.  For  example,  for  second-order 

processes  (0  =2),  n.  =10  if  the  model  tests  are  performed  on  a  scaled- 

r  l. 

i 

down  burner  in  which  all  lengths  have  been  reduced  by  a  factor  of  10.  At 
!-*r  -1  -1 

the  same  time,  n^  *  nf  =  10  ^  np=  ,  i.e.  ,  if  model  tests  have 
been  carried  out  at  densities  ten  times  larger  than  those  which  are  of 
interest  for  the  full-scale  device. 
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Combining  Eqs.  (19)  and  (21)  we  find  that 


kL,  PI 


whence  it  follows  that  simultaneous  similarity  with  respect  to  the  groups 

D,  .  and  r.  can  be  maintained  only  for  first-order  processes  since 

1.1.  r  i  1  * 

we  expect,  in  general,  that  n  oc  n 

El,pi 

possible  to  maintain  simultaneously  similarity  in  model  tests  and  in 
large-scale  burners  with  respect  to  the  important  similarity  groups 

O,  .  and  f  .  except  in  the  unrealistic  case  that  the  effective,  overall 

1.1,  r  i  K 

reaction  order  is  unity. 


To  summarize,  it  is  not 


III  SCALING  OF  FLOW  WITH  MAINTENANCE  OF  CONSTANT  VALUES 

OF  D.  .  AND  Bu  .  FOP  TRANSPARENT  GASES 
x,i, r  *i,i 

The  scaling  procedures  discussed  in  the  preceding  Section  II 
were  designed  to  assure  similarity  of  chemical  reaction  profiles  and  of 
radiative  energy  loss  rate  per  unit  area  relative  to  the  convective  energy 
transport  rate  per  unit  area.  Thus  they  were  designed  to  maintain  the 
flow  field  in  a  reacting  system  invariant  under  the  influence  of  radiant 
energy  loss. 

For  some  applications  it  may  be  more  important  to  require 
invariance  of  the  chemical  composition  profile  and  of  the  absolute  value 
of  the  radiant  energy  emission  rate  per  unit  area  of  reaction  front  since 
the  occurrence  of  radiant  energy  loss  generally  constitutes  only  a  small 
perturbation  on  the  flow  field.  In  this  case,  it  is  pertinent  to  demand 
invariance  of  the  Bouguer  numbers  BUp.  .  as  well  as  of  D.  .  . 
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In  order  to  assure  invariance  of  Bu.,.  . ,  we  impose  the  condition 

i  t  1 


E  t  -  E  T 

L.Pl.o  H  i;H~  L,Pl,o;Mo  i;Mo 


or 


n  n.  =  1. 

el.p>  l‘ 


(22) 


Equations  (18)  and  (21)  apply  as  before  for  n  =1.  If  we  combine  Eqs.  (21) 


and  (22)  we  obtain  the  result 


n  =  n 

EL,  pi  f 


«r-> 


i.e.,  the  condition 


n  =  n 

E  f 

L,P1 


is  now  only  satisfied  for  effective,  overall,  second-order  processes. 

To  summarise,  simultaneous  invariance  of  the  similarity  groups  0T  . 

1 1  T 

and  Bu_  .  is  possible  only  for  the  important  case  of  second-order, 

overall  rate  processes.  For  this  case,  the  groups  Re^  and  M  are  also 

maintained  invariant  for  a  given  combustible  mixture. 

The  foregoing  considerations  lead  to  the  important  conclusion 

that  significant  model  testing  of  the  radiative  properties  of  reacting  gas 

flows  is  possible  for  those  cases  in  which  radiant  energy  emission  does 

not  produce  significant  perturbations  in  the  flow  field  (e.g.,  transparent 

gases)  provided  we  are  satisfied  in  scaling  the  radiative  properties  per 

unit  area  of  reacting  mixture.  The  correlated  changes  in  geometrical 

specifications  (L.)  and  operating  density  (p  )  are  determined  through 

Eq.  (21)  for  n  x  n  s  n  s  1. 

P 
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IV.  SCALING  PARAMETERS  FOR  RADIANT  ENERGY  EMISSION  FROM 
ISOBARIC  AND  ISOTHERMAL  SYSTEMS  FOR  ARBITRARY  OPACITIES 
AND  VARIOUS  SPECTRAL  LINE  AND  MOLECULAR  BAND  MODELS 

For  isothermal  systems  at  the  temperature  T,  the  emitted 

•teradiancy  may  be  calculated  from  the  equation 

B(T)  =  ijR°[l-exp(-PwX)]  du>  . 

The  result*  of  these  calculations  are  summarised  in  Table  la  for  selected 
spectral  line  shapes  and  in  Table  lb  for  selected  vibration-rotation  band 
models. 

From  the  data  listed  in  Tables  la  and  lb,  we  may  draw  the 
following  important  conclusions: 

1.  For  transparent  gases,  the  steradiancy  is  directly  proportional 
to  y  L  irrespective  of  the  spectral  line  contour  or  of  the  band 
model. 

2.  The  important  scaling  parameter  is  f>  L  at  all  optical  depths  for 
isolated  Doppler -broadened  lines,  and  for  all  band  models  in 
which  the  spectral  line  structure  is  effectively  smeared  out. 

3.  The  steradiancy  is  proportional  to  f  ^L  at  moderate  to 
large  optical  depths  (a)  for  isolated,  collision-broadened 
lines,  (b)  for  isolated  lines  with  combined  collision  and 
Doppler  broadening  falling  in  the  "square-root  region"  of 
the  curves  of  growth.  Also  for  statistical  distributions  of 

the  lines  described  under  (a)  and  (b)  "\J  L  is  the  „ 

important  scaling  parameter. 


Table  la.  Radiation  Scaling  Rule*  for  Selected  Spectral  Line  Shapes,  Isothermal  Emitter*. 
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Table  lb.  Radiation  Scalin£  Rule*  for  Selected  Vibration-rotation  Band  Modele,  Isothermal  Emitter 
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risk  indicate  that  the  band  radiancy  is  directly  proportional  to 
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V.  SCALING  PARAMETERS  FOR  RADIANT  ENERGY  EMISSION 
FROM  ISOBAFIC  BUT  NON -ISOTHERMAL  SYSTEMS  FOR 
ARBITRARY  OPACITIES  AND  VARIOUS  SPECTRAL  LINE 
AND  MOLECULAR  BAND  MODELS 


One  of  the  cl&aaical  approaches  to  the  theoretical  calculation 
of  radiant  energy  emission  from  non-isothermal  systems  is  exemplified 
by  the  Lundblad  aeries  development  for  the  solar  photosphere. ^ 

The  spectral  steradiancy  at  the  frequency  v  and  at  the  optical 
depth  T'y  =0,  corresponding  to  the  geometrical  length  s  *  0,  in  the 
direction  01  (see  Fig.  2),  is  given  by  the  relation 


B„ 


oo  . 

(0.0')  =  J  Bv°(rv)|exp[-rvsec0']|  (sec0')dry 


00 

r 

0 

r 

-J  B>)[exp[. 

'J  kL,v<*'H] 

(23) 


where  B^°  is  the  blackbody  steradiancy  for  local  thermodynamic 

s 

equilibrium  at  the  optical  depth  =cos0'J  k  (s')ds'  corresponding 

0 

to  the  geometric  length  s,  along  the  beam  of  the  emitting  system,  for 
a  spectral  linear  absorption  coefficient  ky  =  y  . 

If  B^(  Ty)  is  developed  in  a  (Lundblad)  power  series  in  Ty  , 

Bj’lV-I  W  ■  <«> 

then 

oo 

Bv(0,  O')  =  ^  a.  cos1©'  J  yie’ydy  =  ^  a^'.  cos1©'  (25) 

0  0  0 

where  the  substitution  y=  T  secO'  has  been  used.  Comparison  of  Eqs.  (24) 


i=0 


oo 


oo 
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and  (25)  shows  that 

Bv(0,0')  =  Bv°(Tt^  =co|8,)  for  a.=0  if  i  a  2,  (26) 

i.e.,  if  only  the  first  two  terms  are  used  (Eddington-Barbier  approxi¬ 
mation)  in  the  power  series  given  in  Eq.  (24).  The  physical  interpre¬ 
tation  of  Eq.  (26)  is  the  following:  the  spectral  steradiancy  at  T  =0, 
observed  at  an  angle  O'  ,  for  a  non-i Bother mal  system  is  identically 
equal  to  the  numerical  value  of  the  spectral  blackbody  steradiancy 

Bj*(Tr  _co>qi)  **  the  optical  depth  =cos0'  or  at  the  geometrical 

^  ®  a 

length  defined  by  J  k  ds'  =  1,  provided  only  two  terms  are  used  in 
0  U*Y 

the  power  series  expansion  shown  in  Eq.  (24). 

It  is  interesting  to  consider  the  possible  temperature  profiles 

for  selected  spectral  line  shapes  that  are  consistent  with  the  statement  s 

s 

=  »Q+  ‘iV5  a0+  ‘l00*0’  /  kv(*')d»' 

0 

s 

»  By°(0)+  [Bv°(Tr^=coBQl)-By°(0)]  f  kv(s’)ds'  (27) 

v  Jo 

and 

s 

kv(s')ds*  =  1.  (28) 

We  assume  a  monotone  variation  of  T  and  T  with  s  (and  thus  also  of 
with  T)  and  we  impose  the  boundary  condition  T=Tq  at  =  «  =  0. 
Differentiation  of  Eq.  (27)  with  respect  to  T  yields  the  differential  equation 

Note  that  s  is  measured  along  the  direction  O'  shown  in  Fig.  2. 
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(coa  0')k 


,(T)^ 


_1_ 

a. 


dV 

ST 


_1_ 

a. 


,.2  4 
2h  v 


1  exp(h  w/kT) 

"kT2  [exp(hWkT)-l]4 


(29) 


It  ia  now  possible  to  apecify  kyl(T)  for  varioua  spectral  line  profilea 
belonging  to  varioua  aaaumed  atomic  or  molecular  emittera.  We  may 
then  integrate  Eq.  (29)  in  order  to  find  a  aa  a  function  of  T.  Finally, 

Eq.  (28)  may  be  uaed  to  obtain  the  proper  value  of  a,  and  hence  of  T, 

a 

for  which  B.,(0,  O')  =  B  °(T  for  [  k  (a')da'=l). 

y  0  y 

In  the  analyaia  presented  in  Eqs.  (23)  to  (29),  it  has  been 
assumed  that  the  quantities  a^  are  independent  of  T ^  .  For  an  emitting 
system  with  structure,  this  statement  can  be  true  only  spectrally,  i.e.  , 
a  different  value  of  a^  must  be  chosen  at  a  different  frequency  for  any 
specified  temperature  dependence  on  geometrical  length.  The  implications 
of  this  fact  may  be  clarified  by  referring  to  the  schematic  diagram  shown 
in  Fig.  3.  The  temperature  profile  must,  of  course,  be  independent  of 
frequency  in  any  physically  meaningful  problem.  However,  the  physical 
location  a  and  the  temperature  T  at  which  Eq.  (28)  is  satisfied  is 
strongly  dependent  on  frequency.  In  the  near  line  wing  at  the  frequency 
A  the  integral  condition  of  Eq.  (28)  will  be  met  for  small  values 
of  s  and  T;  on  the  other  hand,  in  the  far  wings  of  spectral  lines  where 
V  =  Vq+ A  v^,  much  larger  values  of  s,  and  hence  of  T,  are  required 
(compare  Fig.  3)  In  other  words,  the  contributions  to  B^fO.O') 
calculated  according  to  Eq.  (26)  arise  from  regions  of  different 
temperature  at  different  locations  for  different  frequencies  in  such  a 
way  that  the  far  line  yrings  will  make  relatively  larger  contributions 
since  they  may  be  "seen"  at  greater  geometrical  depths  and, 
correspondingly,  at  higher  temperatures. 
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Fig.  3.  Schematic  diagram  showing  the  relation  between  tempera¬ 
ture  T  and  diatance  •  for  a  specified  angle  0  .  The  values 
of  s  and  T  required  to  satisfy  Eq.  (28)  are  shown  for  the 
representative  frequencies  vQ  +  Avj  in  the  near  wing  and 
v  +  Av^  in  the  far  wing  of  a  line. 
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Since  &j  may  vary  with  frequency,  it  will  prove  to  be  convenient 
to  introduce  a  frequency  dependence  for  a^  deliberately  in  such  a  way  as 
to  allow  a  universal  representation  of  a  reduced  distance  variable 
(which  is  a  function  of  y  )  as  a  function  of  T. 


A.  Isolated  Spectral  Lines  Belonging  to  Diatomic  Emitters  with  Collision 
Broadening 

For  local  thermodynamic  equilibrium,  we  find  for  diatomic 
emitters,  to  the  harmonic -oscillator  and  rigid-rotator  approximation, 
the  following  relation  for  collision  broadening: 

2  0-  T  ,3/Z 

kv=  Au*l  gu('rf^)B^('^r)  [l-«*(-u0T0/T)][l-exp(-hV0/kT)] 

o  j 

X  [exp(-Ei/kT)]  |l+[(v-Vo)2/(bo2To/T)]|  .  (30) 

Here  vQ  is  the  frequency  at  the  center  of  the  emitted  spectral  line: 

t)  is  the  Einstein  coefficient  for  spontaneous  emission  for  the 
transition  producing  the  given  spectral  line;  gu  is  the  statistical  weight 
of  the  upper  energy  level  involved  in  the  transition;  p/kTQ  represents  the 
number  of  molecules  per  unit  volume  at  the  pressure  p  and  at  the 
reference  temperature  Tq  with  k  denoting  the  Boltsmann  constant; 

0~o=hcB/kTo  where  B  is  the  appropriate  rotational  constant  for  the  rigid 
rotator;  bQ  is  the  spectral  line  semi-half-width  at  the  pressure  p  and 
at  the  reference  temperature  Tq  and  we  have  assumed  that  b=bQ Tq/T 
at  the  constant  pressure  p;  u  =hv  /kT  where  y  represents  the  normal 
vibration  frequency  of  the  diatomic  molecule  (harmonic  oscillator); 


E 


l 


-  energy  of  the  lower  state  above  the  *ero-point  energy. 
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From  Eq».  (29)  and  (30)  we  obtain  the  following  differential 

equation  f  after  approximating  v  in  Eq.  (29)  by  VQj  : 

1 

^  bQ  /  Tq  \I  exp(-hVQ/kT)  exp(E^ /kT) 

***  “  cos©'  a,  (  T  / 


[l-exp(-hVo/kT)]i  [l-exp(-uoTo/T)]  \!+  T>- 

(31) 


(v-Vq)  ' 


where 


,  16.V  "o 


<r  T  g  A 

A  A  ■«  1 


(32) 


1.  The  Limiting  Case  hVQ/kT  <<  1,  u0To/T<<1,  (h  V0*E^  )/kT  « 1 
For  hv  /kT « 1,  u  T  /T«l,  (hv  -E,  )/kT«  1.  Eq.  (31) 

O  O  O  O  4> 

becomes,  after  integration  between  the  limits  8=0  at  T=Tq  and  s,  T, 


s  * 


jA 

cosO' 


(T9/2-T 


9/2 


£ 

coi 


)+£^' 


,,.11/2  _ 

(T  -T 


11/2 


). 


(33) 


where 


A '  =  | 


i> 


- 5 — 77F 

u  (hv  /kT  )  T  1 
o'  O  0  0 


,  2  ^  <bo/al) 

<3  =  Ti 


(34) 


v  -  V 


3  9/2 

u  (hv  /kT)  T 
o'  o  o 


2.  The  Limiting  Case  h  VQ/kT  »  1,  uoTq/T»1 

For  hv  /kT  »  l,  u  T  / T»l,  and  with  x  s  (hv  -E, /kT, 

O  O  O  u  4/ 

Eq.  (31)  becomes  . 


cr 

o 

1 

fkVE/' 

I 

-3/2  -x 

1  + 

2  -E. 

o  i  _-l 

COS& 

r  Sj  o 

VT 

i  °  i 

Z  € 

k°  i 

kT 

o 

dT 


dx.  (35) 
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Integrating  again  from  8=0  at  T =Tq  and  z=zqs  (hVQ-E^)/kTo  to  a, 
we  find  now  that 


(coaQ')a 


A 


1  3 


/  hv  -E . 

I 

v  \T 

V  -  v  \ 

b 

h v  -E. 

I  o  i 

i,+  At  -£■ 

1  o  a. 

o  i  | 

O  j 

\ 

kl 

b 

1 

°  / 

°  / 

x, 


(36) 


where 


z 

h--J 


-3/2  -z  ,  Y 

z  e  dz,  I. 


z 


-5/2  -z. 
z  e  dz. 


(37) 


In  order  to  evaluate  Ij  and  I it  ia  convenient  to  write  the  identity 


z  z  .  o 

/-n  -z ,  f  -n  - z ,  -n  - 

z  e  dz  =  -  |  z  e  dz  +  I  x  e 

Z  00  OD 


zdz  for  n  =  y  or  j,  (38) 


where  the  aecond  integral  appearing  on  the  right-hand  aide  of  Eq.  (38) 
ia  negligibly  amall  compared  to  the  firat  for  zq  much  larger  than  z.  But 


md 

•/ 

OD 


-3/2  -z, 
z  c  dz 


2  e-  ,  r(  1  r*(2-M 

.2_.2r(r)h. 


(39) 


where  P  ( j)  ia  the  complete  gamma  function  of  argument  j  and  rj  y) 

1  8 

ia  the  incomplete  gamma  function  of  z  of  argument  y  •  Making  uae  of 

g 

the  notation  of  Pearaon,  the  preceding  expreaaion  may  be  rewritten  in 
the  form 


z 

-/ 

oo 


-3/2  -z. 
z  e  dz  =  2  i  ■ 


r  e  ‘ 


(  vr 


-  VT[1-I(. 


VoTs 


■0.5 


!)] 


(40) 


With  the  available  table  a  of  the  incomplete  P  -function,  which  give  a 
valuea  to  aeven  aignificant  figurea,  Eq.  (40)  can  only  be  evaluated  for 


8  K.  Pearaon,  Tablea  of  the  Incomplete  T  -Function,  Cambridge 
Univereity  Preaa,  Cambridge,  195^. 
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x<9.  For*  >9,  the  integral*  may  be  evaluated  either  numerically  or 
elee  by  ueing  a  aimple  approximation  procedure. 

Integrating  by  parte  twice  yields  the  expreeeion 
z  z  z  z 


-n  -z  , 
z  e  dz 


-n  -* 
z  e 


'+  j  J  z  ne  Zdz  +  J  ( J  *  n***d*^)d*‘  (41) 


00 


00 


00  00 


For  sufficiently  large  values  of  z,  Eq.  (41)  reduces  to 


-n  -z 

-  -p+OTTJ 

00 


■ 


-n  -z  , 
z  e  dz 


(42) 


or,  in  somewhat  cruder  approximation, 
z 

/*n  j  * n  ~ z  . . > . 

z  e  dz  x  z  e  (43) 

ao 

Results  obtained  by  using  Eqs.  (42)  and  (43),  and  also  by  using 

g 

Eq.  (40)  together  with  tabulated  values  of  the  incomplete  gamma  function, 

are  plotted  in  Figs.  4a  and  4b.  These  data,  together  with  Eqs.  (36),  (37) 

and  (38)  yield  the  desired  temperature  profile.  For  hv^/kT  «  1,  the 

dependence  of  T  on  s  is  easily  computed  by  using  Eq.  (33). 

Reference  to  Eqs.  (33)  and  (36)  shows  that  the  temperature  profile 

depends  on  the  frequency.  At  the  line  center,  however,  (8  =0.  Therefore 

s  may  be  computed  as  a  universal  function  of  T  for  hvo/kTQ<<  1. 

Similarly,  for  h  v1  /kT  >>  1,  the  second  term  in  Eq.  (36)  vanishes  and 
o  o 

(cosO')  s/|i4(bQ/a1)To^(hvo-E^  )/kTQj^2j  =  Ij,  which  has  been  plotted 

in  Figs.  4a  and  4b  for  the  special  cases  zq=  ao  and  zq=  20:  in  Figs.  5a  to 

5c,  the  corresponding  temperature  profiles  are  shown  for  zq=  20, 

T  =  300°K. 
o 
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I*  a  7.'1 

For  comparison  also  - '  -  dz  is  plotted,  which 

-00  (1-e-V 

determines  s(T)  for  z  pf  the  order  of  1  [see  Eq.  (31)]. 


-2 


10  12  14  16  IS  20 


e*z  z'ndz  as  a  function  of  z  for 
oo  - 

10  <  z  s  20  (z  =  20). 
o 


A  plot  of  the  firat  approximation  to  1^  [aee  Eq.  (43) J  aa  a  function 
of  T  and  a  (for  zQ  =  20,  J0  *  300°K,  E^  *  0).  The  contribution 
of  the  term  e  a  "  tq  I.  ia  negligibly  amall. 
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In  the  line  wings,  the  firat  terms  of  Eqa.  (33)  and  (36)  become 
negligibly  email.  Hence  reduced  temperature  profiles  are  again  determined 
in  terms  of  easily  computed  quantities  or  in  terms  of  I^. 


3.  Determination  of  Steradiancy  B(0,  0)  for  hVQ/kTo  »1,  uqTo/T>>1 

If  a  temperature  profile  is  specified  which  differs  from  Ij  or 
by  at  most  a  constant  factor,  then  the  steradiancy  B(0,0)  may  be  calculated 
for  suitable  frequency  regions.  From  Eq.  (26)  it  is  apparent  that 

oo 

B<0,0)  =  J  BV°(T  =co,0,)dv  ,  (44) 

0 


where  Bv°(Tr^ 
Tv/cosQ'  =  1. 


_cOBgi)  the  blackbody  steradiancy  at  that  location  where 
According  to  Eq.  (27), 


al  = 


V(TTv  =co.e,)  -»y  (V 


/  cosO' 


(45) 


Let  us  consider  the  temperature  profile  (cosO1)  s=LjIj,  where  the 
characteristic  length  Lj  is  a  constant.  Close  to  the  line  center,  Eqs. 
(36)  and  (45)  lead  to  the  relation 


V*Tr„  *co,0')  -Bv(To)=  I"  bc 


TT 


T  cosO1 
o 


(46) 


Next  we  introduce  Eq.  (46)  into  Eq.  (44)  and  integrate  from  V  -  A  v>c  to 

V  +Av  ,  where  A  V  is  sufficiently  small  to  justify  use  of  the  first  term 
o  c  c 

only  in  Eq.  (36).  In  this  manner  we  obtain 
v  +av 


/ 


<4t 


B>rw=Co.0-^“  2Tf  bd 


1  hvo-E£ 
— JcT 

\ 


A  cosO',  (47 ) 
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where  Bj^T  )  has  been  neglected.*  The  contribution  of  the  term  in 
Eq.  (38),  which  has  been  neglected  in  Eq.  (47),  is  close  to  (A  vo/bQ)^. 

Let  us  now  consider  a  temperature  profile  (cosO1) 
and  investigate  only  the  contribution  of  the  wings  in  Eq.  (36).  Then, 
proceeding  as  before, 


'hv  -E  \ 
o _ I 


kT 


2!  v'-v' 


°]  cosO*  +  B^(T  ).  (48) 
'o  / 


Introducing  Eq.  (48)  into  Eq.  (44),  we  obtain  the  following  integral 
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o  w 
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(49) 


where  By(TQ)  has  been  neglected.  The  integration  limit  towards  the  line 
center.  vQ+  must  be  chosen  sufficiently  large  so  that  the  first  term 

in  Eq.  (36)  is  negligibly  small.  The  upper  limit,  w  must  be 


consistent  with  the  restriction  that  Ty  =cosO'  for  all  v-1  .  Clearly  the 

basic  relation  given  in  Eq.  (26)  can  only  hold  provided  T_  £  T  , 

Ty  =cosw  max 

where  T  is  the  highest  temperature  of  the  system,  i.e.,  the  value  of 
max 

is  determined  by  the  expression 


AV  u 
w 

’  Bv°<T„,„>4 

kT.  1 

M (  W 

- - 1 

bo  ' 

K -E!  j 
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j 

(50) 
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B^°(Tq)  may  be  neglected  ordinarily  unless  T  -  constant  or  cosO'  -  0. 
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B.  I •olated  Spectral  Lines  Belonging  to  Diatomic  Emitters  with 
Doppler  Broadening 

For  Doppler -broadened  lines,  we  find 


k.,= 


8ir  v 


1  *u-.X  gu  Tt" 


Xn2 

I 

fTol 

I 

r 

bD,o 

n 

[i-.«P(-u0T0/t] 


X  [  l-exp{-hvo/kT)][  exp(-E^  /kT)]  exp 


V-V  ,  T 


(54) 


where  the  Doppler  half  width  under  reference  conditions  is  given  by 

/  2kV”2  \ 1/2 


’d,o - r 

me 


From  Eqs.  (29)  and  (54)  we  find  that 
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where  i*  given  by  Eq.  (32). 


1.  The  Limiting  Case  hv^/kT  <<  1,  uqTo/T  <<  1,  (h  v^-E^  )/kT  <<  1 


o  “i 


Equation  (55)  becomes  now 

1  .  .9 
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We  define  the  region  near  the  line  center  by  the  condition  that 
the  exponential  term  may  be  replaced  by  unity.  In  this  case 
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(cosO')s  =  Jc  (T/Tq)11/2-1 


where 


i  n  .  JL  D»P  0 

**■  "11  u  a. 

O  1 


* 

ir£n2  I  hv  i 


When  the  exponent  ia  sufficiently  large,  the  following  approxi¬ 


mation  may  be  used: 


(cosO-')s  =  £  1  exp  (InZ) 
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2.  The  Limiting  Case  hv  /kT  »  1,  u  T  /T  »  1 

o  o  o 

Integration  of  Eq.  (55)  yields  the  expression 
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For  w  >  0,  we  may  proceed  as  with  the  collision-broadened  line  and 

the  integral  w 

f  -5/2  -w.  .... 

I,  =  -  /  w  e  dw  (60) 

w 

o 

may  be  evaluated  by  using  the  methods  of  Section  V,  A(2). 

Having  assumed  hvQ/kTQ>>  1,  the  frequency  region  for 
w  <  0  occurs  far  out  in  the  line  wings  and  contributes  relatively  little 
to  the  total  steradiancy. 


3.  Determination  of  Steradiancy  B(0,  0)  for  hv^/kT  >>  1,  uoTq/T»  1, 
(hVQ-Ex)/kT  »  1 


The  frequency  dependence  of  ds/dT  given  in  Eq.  (59)  does  not 

permit  us  to  choose  a^V)  in  such  a  way  that  a  reduced  distance -temperature 

profile  can  be  constructed.  The  difficulty  is  caused  by  the  occurrence  of 

a  product  of  frequency-  and  temperature -dependent  terms  in  the  exponent. 

Near  the  line  center,  we  may,  however,  calculate  the  radiant 

flux  since  the  integral  1 2'  is  determined  almost  entirely  by  its  upper 

limit  w  (corresponding  to  T  =  T  )  provided  that  w  -w  >2.  This 

min  max  o  min 

property  of  the  integral  has  been  discussed  in  Section  IV,  A(2). 

Introducing  the  additional  restrictions 
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Eq.  (59)  becomes 
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where 


I. 


Z 
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ao 


z  ^2  e  Zdz,  z  =  (hv  -E.  )/kT, 

o  X 


and  z  ia  independent  of  V  . 

The  inequalities  in  Eq.  (61)  become,  for  typical  fundamental 

vibration-rotation  bands  of  diatomic  molecules, 

2 

«  30 


and 


fv-*'oV 

/  v-v  \2  T 

o  i  „  ,  ,  max 

*TT~  j  «  1  • 4  ^  . 

D,  O  ;  o 

Hence,  for  large  values  of  T _ /T  ,  the  correct  integral  is  obtained  for 

max  o 

a  frequency  range  that  may  be  appreciably  larger  than  b_  .  On  the 

Uf  o 

other  hand,  for  systems  with  small  temperature  gradients,  we  have 
obtained  a  temperature  profile  that  is  applicable  only  very  close  to  the 
line  center. 

For  the  temperature  profile  s(cosQ')=  L^I^,  we  find,  after 
integrating  over  the  frequency  range  V0-Av^to  £  compare 

Eq.  (47)]  . 


that 


/ 


*Vavd 


Bv<Trv=co.0.>d^2 


■^b_  T  ,  , 

D,o  o  1 


TnT 


^/hy 

kT 


(cosO') 


AV 


D’ 


(63) 


-38- 


In  general,  the  integral  represents  only  the  contributions  arising  from 
a  narrow  frequency  range  near  the  line  center.  Since  it  is  reasonable 
to  assume  that  AV  is  proportional  to  b_  ,  the  radiancy  for  this 
frequency  range  near  the  line  center  is  inversely  proportional  to 


C.  Gray  Body 
With 


k_  (T)  =  k  =  constant, 


integration  of  Eq.  (29)  leads  to  the  expression 


s(cosO') 


1  2hV 
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a^C  c 


1 _  1 

exp(hv/kT)-l  ~  exp(hv/kTQ)-l 


(64) 


In  certain  regions  of  frequency  and  temperature,  Eq.  (64) 
reduces  to  a  universal  relation  between  s  and  T. 


1.  The  Special  Case  hWkT  <  hw/kTQ  «  1 


Equation  (64)  now  reduces  to 

s(cosO')  =  -~^(T-T  ). 
a.kc 


Hence,  for 

Eqs.  (44)  and  (45)  lead  to 
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s(cosO')  =  L  (T-To), 


(65) 


(66) 


J  By(0,Q)<h>  =  J 
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By°(To)dV+  -y  (V23-V13)(cosQ').  (67) 

3  Lite 


r  c  —  -l 

and,  if  T  »  Tq,  then  J  Bv(0,  8)dvoe  (  o  L  )  ;  the  scaling  parameter 

V1 

is  p  L  . 

Jo 
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2.  The  Special  Case  hWkTQ>>  1,  hv/kT«  1 

In  Eq.  (65),  T  replaces  (T -T  )  and  the  preceding  reaulta  apply. 

3.  The  Special  Case  hv/kTQ  >>  1,  hv/kT  >>  1 

In  this  case  there  is  no  s-T  curve  which  is  independent  of  v  . 

VI.  RADIATIVE  SCALING  PROPERTIES  FOR  REPRESENTATIVE 
TEMPERATURE  PROFILES 

The  integral  expression  of  Sec.  I.  C.  for  the  spectral 

steradiancy  B v  is  the  formal  solution  to  the  linear,  first-order 

differential  equation 

dBv 

- w  =  Lokv(Bv°-Bv  )  (68) 

ds 

where  s*  =  s/Lq,  is  the  distance  along  the  line  of  sight,  and  LQ  is  a 
characteristic  length  of  the  system.  Equation  (68)  has  been  integrated 
numerically  by  means  of  a  fourth  order  Runga-Kutta  method  for 
representative  temperature  profiles.  The  temperature  profiles  are 
represented  by  the  expressions 

T  s  (Tmax'To)(l*l,*"l|m,+  To*  mal*2'4-  <69> 

The  specified  temperature  profiles  are  sketched  in  Fig.  6. 

The  absorption  coefficient  k^  is  given  by  Eqs.  (30)  and  (54) 
for  dispersion  and  Doppler -broadened  line  contours,  respectively.  The 
spectral  and  total  line  steradiancies  at  s  *  2  have  been  computed  for  a 
typical  strong  line  of  the  vibration-rotation  spectrum  of  the  hydrogen 
fluoride  molecule.  We  have  chosen  the  values  Tq  =  300  K  and 
T  ■  3000°K.  Representative  results  of  the  calculation  are  shown 
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in  Figa.  7,  8,  and  9  for  the  R3(v»0,  J«3  "•  v«l,  J*4)  lino  of  HF.  For  this  line 

in  case  of  dioperaion  broadening.  bQ=  .  132p  cm”'  (p  in  atmoe)  whereaa  for 

-3  -1 

Doppler ‘broadening  we  have  uaed  bQ=  5.  85  x  10  cm  The  reference  optical 

depth  ia  defined  aa  t  =  (L  /wb  )f  k.  dv  (k  =abaorption  coefficient  evaluated 
o  o  o  *  v»o  v»o 

at  T0).  Bv  ia  the  apectral  ateradiancy  at  a*«  2.  and  B>#  repreaenta  the 
integrated  ateradiancy  at  a**  2  for  the  entire  line. 

Reference  to  Fig.  9  ahowa  that  Doppler  and  diaperaion  broadening 
produce  the  aame  total  ateradiancy  for  the  caae  under  conaideration  for 
large  value,  of  provided  that  bodi<p#r#ion  -  2b0>  DoFpl„.  i.  ••  ,  for 
p  **  0.  1  atm. 

A.  Diaperaion-Broadened  Line 

Examination  of  Fig.  7  ahowa  that  the  R3  line  of  HF  retaina  a  typical 
diaperaion  contour  until  becomea  greater  than  about  10.  when  eelf-reveraal 
become,  important.  Hence  ro  “■  10  may  be  a  aid  to  define  the  upper  limit  of  the 
traneparent  gaa  regime.  From  Fig.  9  it  ia  aeen  that  the  preaaure  and  length 
dependence  of  the  ateradiancy  for  to  <  10  ia 

».*  *  P0Lo'  (70) 

1.  e.  ,  it  ia  the  aame  aa  for  an  iaothermal  traneparent  gaa.  In  general,  tq*  10 

correaponda  to  a  email  phyeical  length  for  a  atrong  apectral  line.  For  the  R3 

line  of  HF,  the  value  of  L  at  t  =  10  ia  about  0.  015  cm. 

o  o 

The  center  of  the  line  ia  eaaentially  completely  aelf-abaorbed  for 

t  “■  100,  *  Figure  9  ahowa  that  for  r  >100 
o  o 

*The  condition  that  B  #  «(p0^L0)^  for  t  exceeding  the  value  required  to  make 
the  line  center  "black”'  may  be  uaed  to  derive  an  approximate  relation  for  the 
critical  minimum  value  of  t  above  which  Eq.  (71)  appliee  for  varioua  tempera¬ 
ture  profilee  [for  detail.,  ale  the  Ph.  D.  theaie  of  M.  Thomaa,  California 
Inetitute  of  Technology,  Paeadena,  California,  June  1964]. 
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Fig.  7.  The  spectral  steradiancies  for  m  =  2  for  the  R3  line  of  HF  at  s  = 

as  a  function  of  (v  -  vo)/bQfor  dispersion  broadening  and  various  values 
of  the  reference  optical  depth  T  . 


(I0-14  watt /cm*) 


Fig.  8.  The  spectral  steradiancies  B  *  for  m  =  2  for  the  R3  line 

♦  v,  s 

of  HF  at  s  =  2  as  a  function  of  (v  -  vQ)/bo  for  pure  Doppler 
broadening  and  various  values  of  the  reference  optical 
depth  T  . 
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B.*0C  V^o2C  '  <7l> 

which  ia  the  aame  as  for  the  Eddington-Barbier  approximation  or  for 
the  iaothermal  caae  for  large  optical  deptha.  For  a  atrong  line,  the 
tranaition  to  the  regime  deacribed  by  Eq.  (71)  occura  at  a  amall  phyaical 
length  (  —  0.15  cm  for  R3  line  of  HF) 

B.  Doppler -Broadened  Line 

The  plot  in  Fig.  8  ahowa  that  aelf-reveraal  for  the  Doppler- 
broadened  line  becomea  important  for  as  5.  However,  reference  to 
Fig.  9  indicatea  that  Eq.  (70)  remaina  valid  up  to  To  s  10  which,  at 
psO.l  atmoa  ,  correaponda  to  LQ  ^  0.007  cm  for  the  R3  line  of  HF.  For 
larger  valuea  of  T  ,  however,  the  Doppler -broadened  line  doea  not 
approximate  the  behavior  of  the  diaper aion-broadened  line,  i.e. ,  B  does 
not  become  aimply  proportional  to  a  power  of  Lq.  Rather,  Bbi>  become*  a 
weaker  and  weaker  function  of  T  aa  aaturation  ia  approached. 

The  more  complicated  behavior  of  the  Doppler -broadened  line 
compared  to  the  diaperaion  broadened  line  ia  the  reault  of  the  fact  that 
the  Doppler  line  half-width  increaaea  with  temperature  while  the  diaperaion 
line  half-width  decreaaea  with  temperature  at  conatant  preaaure.  Hence, 
viewing  a  Doppler -broadened  ayatem  at  a  =2 ,  the  radiation  emitted  from 
the  higher  temperature  regiona  ia  "Been"  at  all  optical  deptha  aince  thia 
radiation  ia  not  extenaively  reabaorbed. 


